An algebraic approach to problems with polynomial Hamiltonians on Euclidean spaces 
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Explicit expressions are given for the actions and radial matrix elements of basic radial observables 
on multi-dimensional spaces in a continuous sequence of orthonormal bases for unitary SU(1,1) 
irreps. Explicit expressions are also given for SO(7V)-reduced matrix elements of basic orbital 
observables. These developments make it possible to determine the matrix elements of polynomial 
and a other Hamiltonians analytically, to within SO(iV) Clebsch-Gordan coefficients, and to select 
an optimal basis for a particular problem such that the expansion of eigenfunctions is most rapidly 
convergent. 



INTRODUCTION 



in 
o 
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This paper is concerned with the large class of problems lying between the relatively few that are exactly solvable 
O , and others that are only solvable by numerical methods. In particular, algebraic methods are developed for computing 
' the spectral properties of Hamiltonians which arc polynomials in the Cartesian position and momentum {pi 

>n 

■ are also considered. Restriction to this sub-class of Hamiltonians is because, for them it is appropriate to employ 
bases of harmonic oscillator or modified oscillator wave functions. Probably a similar treatment can be developed for 
hydrogenic and modified hyrdrogenic bases which will enable the admission of potentials, like the \/r potential of the 
^ \ hydrogen atom, which are not rational functions of the Cartesian coordinates. 

The most widely studied Hamiltonians of this kind are for exactly solvable central force problems. Some central 
force problems are solved by use of various realizations of an SU(1,1) spectrum generating algebra (cf., for example, 
Refs. 0,0,01 for reviews). Others are solved by the so-called factorization method introduced by Schrodinger @ 
and employed extensively by Infeld and Hull [6j • (A pedagogical review of the factorization method is given in several 
\ texts on quantum mechanics, e.g. Q.) Many special case solutions for power law and inverse power law potentials 
have also been obtained by the methods of Refs. d, El El EI- I n an attempt to identify the underlying symmetries 
of exactly solvable systems, Gendenshtein fl2T | introduced the criterion of shape invariance and a connection with 
O j. supersymmetry. Exactly solvable central force problems play an essential role in the current investigation. However, 
their purpose will be to provide suitable bases in terms of which the matrix elements of a much larger class of 
Hamiltonians of interest can be determined algebraically. 

Underlying the solvability of a central force problem is the factorization of its Hilbert space into a product of radial 
and orbital subspaces 



£ 2 (R N ) ~ £ 2 {R+,dr) <g>£ 2 (SV-i), (1) 
where R + is the positive half of the real line, 

£ 2 (R+,dr) = |i?:r->C J \R(r)\ 2 dr < ooj , (2) 

and the elements of C 2 {Sm-i) are square intcgrablc functions on the N— 1 sphere, Sn-i — SO(N)/SO(N — 1), relative 
to the SO(-/V)-invariant measure. Thus, because the Hamiltonian for a central force problem is SO(iV)-invariant its 
Schrodinger equation reduces to a one-dimensional radial equation. 

More generally the radial R + and orbital Sn-i dynamics are coupled by centrifugal forces. Nevertheless, the 
SU(1, 1) x O(N) dynamical group associated with the above factorization continues to provide powerful algebraic 
(and hence exact) methods for determining the matrix elements of polynomial Hamiltonians in suitable bases. This 
is of enormous practical value even if it is ultimately necessary to resort to numerical methods for diagonalizing the 
Hamiltonian matrix. Moreover, it is a huge advantage if bases functions can be found which are close to the desired 
solutions so that the dimensions of the matrices that have to be diagonalizcd to achieve a given level of accuracy are 
relatively small. 

With such application in mind, Armstrong [l^, Haskell and Wybourne [l4|, for example, derived expressions for 
some radial matrix elements in a basis of eigenstates of the three-dimensional harmonic oscillator. It turns out that, 
for many purposes, other matrix elements are needed and are also needed for higher-dimensional spaces. For example, 
the nuclear collective model is defined on a five-dimensional space ^1 El El . And, while harmonic oscillator bases 
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are appropriate for spherical vibrational systems, they are not the most appropriate for rotational- vibrational systems 
such as diatomic molecules and non-spherical nuclei. 

The current paper gives analytical expressions for the matrix elements of needed observablcs in a continuous 
sequence of orthonormal bases for unitary SU(1, 1) x 0(A) irreps. These irreps include, but are not restricted to, 
the irreps of the harmonic (discrete) series. Nor are they restricted to a three-dimensional space. Matrix elements of 
observables that lie in the SU(1,1) Lie algebra are obtained by algebraic methods. Factorization methods are used to 
compute the actions and matrix elements of other needed observablcs. These developments make it possible to select 
an optimal basis for a particular problem and continue to retain analytical expressions for matrix elements as provided 
by the harmonic oscillator bases. Basis functions for £ 2 (Sn-i), which span 0(A) irreps, are commonly referred to as 
SO(A)-spherical harmonics. They are known for SO(3), SO(4) ~ SU(2) x SU(2), SO(5) [I101, and SO(6). However, 
as the following shows, explicit expressions for the wave functions are not needed for the computation of matrix 
elements. What more is needed, beyond the expressions given in this paper, are the appropriate Clebsch-Gordan 
coefficients. These are now available in SO(3)-coupled bases for all A < 6 (cf. section ??). 



II. REPRESENTATIONS OF SU(1, 1) x SO(7V) 



The SU(1,1) Lie algebra (more precisely its complex extension) is spanned by operators {S±, So} which satisfy the 
standard commutation relations 

= 2S , [§o,S±]=±§±. (3) 
A basis for an SU(1,1) irrcp is given by an orthonormal set of states {|A^); v = 0, 1, 2, . . .} which satisfy the identities 

So\\v) = £(A + 2i/)|Ai/>, 



S+\Xv) = + v)(v + l)\\,v + l) 

S-\Xv) = yJ{\ + V-\)v\\,V-\). 



For such an irrep the Casimir invariant 
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Csuii — (So) — —(S+S- + S-S+) — So(So — 1) — S+S- 



(4) 



(5) 



takes the value A(A — 2)/4, i.e., 



Csmi\Xu) = -X(X-2)\Xu). 



(0) 



The SO (AT) Lie algebra is spanned by angular momentum operators which are antisymmetric, £.y = —Cji, 

and satisfy the commutation relations 



[Cij , Cki] — — i 



SjkCu — SuCkj 



The SO (A) Casimir operator is defined by 



(7) 



(8) 



><j 



If an SO(A) irrep, labelled by v, has basis states {\vm)}, then a basis for an SU(1, 1) x SO(A) irrep (A, v) is given 
by the product states 



\Xv;vm) = \Xv) x \vm). 



(9) 



A. Harmonic series representations 



An explicit realization of the SU(1,1) Lie algebra is given in terms of harmonic oscillator raising and lowering 
operators 
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by the 0(-/V)-invariant operators 

where c • c = c,Cj and n = J^cjcj is a harmonic oscillator number operator. Similarly the SO(iV) angular 
momentum operators take the form 

Cn = -i(cJ Cj - ctcj) . (12) 

When the SU(1,1) and S0(7V) operators are realized in this way, their Casimir invariants are given, respectively, 
by 

4Csuii = h(n + N - 2) + ^N(N - 4) - (c f • c f )(c • c) (13) 

and 

\ = n(n + N -2)-(c f -c^^-c). (14) 

It is seen that when acting on the Hilbcrt space C 2 (R N ), according to the above-defined realizations, these Casimir 
invariants are linearly related by the identity 

A = 4(7 S uii-|iV(iV-4), (15) 

in accordance with a well-known duality relationship |2fJ | (also called complementarity |2l|V This duality implies that 
the labels A and v are in one-to-one correspondence. Thus, if the state \XO;va) satisfies the equations 

§-\\0;va)=0, n\XO;va) =v\XO;va), (16) 

then A takes the value 

X = v + N/2, v = 0, 1, 2, ••• (17) 

Moreover, the eigenvalue of A 

A\XO;vm) = v(v + N-2)\\0;vm) (18) 
signifies that v, often referred as seniority, is an SO(iV) angular momentum quantum number. 



B. Modified oscillator representations 

The discrete harmonic series of SU(1,1) irreps which have duality relationships with SO(A^) irreps are invaluable 
for many purposes. However, there is a continuous series of irreps, which are often more useful. These more general 
realization of the SU(1,1) algebra are obtained by expanding the above operators, using eqns. (|10f> . 



S± 



2r ^ + N 



So 



[-V 2 + r 2 ], 



(19) 



where r 2 = ^ x 2 , and V 2 = J2i d 2 /dx 2 is the Laplacian on C 



2/ntJVl 



The Laplacian has the well-known expansion 



(20) 



where 



1 d 

r.N-1 ^ 



,N 



dr 



(21) 
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In this form the SU(1,1) operators act as differential operators on the wave functions for the basis states {|Az^; vm)} 
of equation These wave functions are conveniently expanded 



r.N-1 



(22) 



dr duj is absorbed into the 



*Xv V m(T) = r-( N - l V 2 KZ(r)y vm (u;), 

where, for convenience, the factor of r*" 1 in the volume element dv{r,Lo) 
definition of the wave functions. Thus, {TZ*; v = 0, 1, 2, . . .} is a basis of radial wave function for £ 2 (R + , dr) and 
{3^m(<^)} is a basis of spherical harmonics for £ 2 (Sn-i) defined as eigenfunctions of A 

Ay vm =v(v + N -2)y vm . (23) 

With the Laplacian in the form V 2 = A — A/r 2 , the action of the SU(1,1) operators of equation on the wave 
functions of eqn. (|22|l gives 



with A = v + N/2 and 



( r> 0) =r -(iV-l)/2 



4 A) ^ A 



(r)y m n(u), 



s 



(A) 

± 



(A-3/2)(A-l/2) 



dr 2 



r 2 =F 2r 



(A-3/2)(A-l/2) 



c?r 2 



(24) 
(25) 

(26) 
(27) 



Although the S^ x ' operators have been derived for the discrete values of A = v + N/2, they nevertheless satisfy the 
SU(1,1) commutation relations 



[<si A) ,<sj A) ] 



25, 



(A) 



$ A) ,4 A) ] 



s 



(A) 



(28) 



for any value of A. Moreover, for any real A > 0, they have a unitary representation and define an orthonormal basis 
{11*} for the Hilbert space C 2 (R+) such that 



V(A + !/)(!/ + 



(29) 



For each A > 0, the operators span a specific SU(1,1) representation. This is evidenced by the fact that, for each 
realization, the SU(1,1) Casimir operator acquires the unique value A(A — 2)/4, i.e., 



4un = " 1) " 4 A) ^i A) = |A(A - 2)1, 



(30) 



where / is the identity. 



III. THE RADIAL WAVE FUNCTIONS 

The wave functions for arbitrary A > are derived as for the standard harmonic oscillator and are given by 



Tlt(r) = i-iyJf™-^ L (A-D (r2) e -rV 2) „ = o, i, 2> . . . (31) 

The lowest-weight functions {TZq} are plotted in Fig.Q] 

An important characteristic of the radial Hilbert space C 2 (M + , dr) is that it is independent of the dimensionality of 
the Euclidean space in which it is embedded. Thus, when 2A is an integer, these functions are equal to the radial wave 



75 12 3 4 5 6 7 
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FIG. 1: Radial wave functions {IZq} plotted as functions of r for a range of values of A. (Computed by P.S. Turner [T^-l 

functions for various TV-dimensional harmonic oscillators. For example, the wave functions of the simple harmonic 
oscillator are given by 

u 2n (x) = -±=K /2 (x) = (-1)" = L(-W(x 2 )e-» 2 / 2 , oo < x < oo, (32) 

V2 Vv 7r ( 2n )' 

1 on+l/2„| 

U2)l+1 (x) = -^# 2 (s) = (-1)" 7= 1 = 4 l/2) (^ 2 )^ r/2 . oo < x < oo. (33) 
Radial wave functions for the spherical (N = 3) harmonic oscillator are given by 



u nl (r) =K+*/\r) = r(n *f + 3/2) 4' +V2) (r 2 ) e^ 2 . (34) 

In general, the radial wave functions for an TV-dimensional harmonic oscillator of SO(iV) angular momentum v are 
given by 

u nv (r)=ni +N / 2 (r). (35) 

These identities reflect the fact that harmonic oscillator radial wave functions belong to SU(1,1) irreps which are 
determined by the angular momentum of the accompanying spherical wave function. For example, eigenstates of the 
three-dimensional harmonic oscillator of angular momentum / have wave functions which in spherical polar coordinates 
take the form 

^ n i m {r 7 7 ip) = -u„t(r)Yi m (6,ip), (36) 
r 

where Yi m is a spherical harmonic. 

When 2A is not an integer, the basis functions {TL^} are not the radial wave functions of any harmonic oscillator. 
Nevertheless, they have algebraic properties that are almost as simple as those which are. For many systems, especially 
models that exhibit rotational as well as vibrational states, it is appropriate to use {Rq} radial wave functions with 
A > v + N/2. Such wave functions are eigenfunctions of a modified oscillator Hamiltonian introduced by Davidson 
[22^ for the purpose of describing the rotational-vibrational states of a diatomic molecule. In this model, the atoms 
of the molecule are 'pushed apart' by enhancing the centrifugal potential by the substitution 

W + i) Ki + + 4 r _ 7 v 
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this is equivalent to replacing the harmonic oscillator potential r 2 by a potential 

V(r) = 4 + r 2 (38) 

which has a minimum at r — tq. A five-dimensional version of the Davidson oscillator has been considered as a model 
for nuclear rotations and vibrations [2^, l24l l25j . 



IV. RADIAL MATRIX ELEMENTS 
A. Matrix elements obtained from irreps of the SU(1,1) Lie algebra 
Claim 1: For A > 1, 



r 2 nl(r) = ^{X + v-iyn^r) + \/(A + v ){y + 1) K x u+1 {r) 



(39) 



A-l V/"!r(A + ^ > l j 



/j.>f 



(-iy— m!t(a + z,)^ a 



A - 1 V i/! T(A + /i) 



dr 



-(A + 2v)K$(r) + (A - 3/2)(A - l/2)-^(r) . 



v 2 7zi(r)y„ m (w)=y OTl H 



^) + ^-2^) + (A - 1)2 - ( r^" 1)2 



(40) 

(41) 
(42) 



The first of these equations follows from equation (|29|l and the observation that r 2 is an element of the su(l,l) Lie 
algebra with expansion 



S^ + S W 



25, 



(A) 



Equation 1401) is obtained by rexprcssing equation (|39|l as the recursion relation 



(A + 2v) /* + y/{\ + v + V(* + »)(v + l)fL+i = <W , 



where 



Starting with the value of = 1/(A — 1) obtained from equation 1)72(1 . this equation has solution 



f X = 



(-iy- u ulT{X + ii) 



for n < v . 



A-i y /i! r(A + v) 

Equation l|41|) follows from the identity 

d 2 (A _| )(A _1) _ 

The expression for V 2 is obtained by recalling that, from equations l(19|l . (|24|l and (|25(l . 



25, 



(A) 



(43) 
(44) 
(45) 

(46) 

(47) 



(A/x; i?| V 2 1 Af; u) = (A/i 
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25, 







5| A) + 5 (A) - 25< A) + (A-^(A-l)-(^+^-^)(^+^|) 



\Xv) 



(48) 



Matrix elements of higher even powers of r are given by repeated use of equation l|39|) . 
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B. Matrix elements obtained by the factorization method 

Similar equations for the operators r, 1/r, and d/dr, are obtained by the factorization method Let A(X) 

and A^(X), where X is a real number, denote the operators 

A(X) = ^- + -+r, At(X) = ~ + -+r. (49) 
dr r dr r 



Then 



It follows that 



A(X)At(X) = -j T + X(X ? 1} +r 2 +2X + L (50) 
dr z r z 

A(XM(X) = ~ + X(X 2 + 1] +r 2 + 2X- 1. (51) 



A(X- |)At(A- |) = 4<^ A) +2A, (52) 

At (A- ±)A(A- i) = 4S< A+1) + 2A-2, (53) 

A(-A + |)At(-A+|) = 45^ A) -2A + 4, (54) 

At(-A + f)A(-A + §) = 45 { ^ A " 1) - 2A + 2. (55) 



Equations (|q"2|) and lt55|) give 



A(A- f)At(A- §)ft A = 4(A + i/)ft A , (56) 
At (A - |)A(A - i)^ +1 = 4(A + (57) 



and hence 



At(A - |)A(A - |) [At (A - 1)7^] = 4(A + v) [At (A - £)7#] , (58) 

A(A - |)At(A - i) [A(A - i)^ 1 ] = 4(A + v) [A(A - . (59) 

From these and similar equations obtained from equations 154(1 and (|55|l . it follows (with a choice of relative phase) 
that 

At (A- \)n x v = 2VaT^ a+1 , (60) 

A(A-i)^ +1 = 2VAT^7^, (61) 

At (—A + l)K$ = 2V^TT Tll-\, for A > 1 (62) 

A(-A + l)Klz\ = 2V^TT Tli, for A > 1. (63) 



Thus, it is seen that the A and At operator act as raising and lowering operators for the radial wave functions in close 
parallel with the way the standard raising and lowering operators act on harmonic oscillator states; the parallel is 
exhibited in fig. [21 This observation is exploited in Section IVTl to obtain the actions of the harmonic oscillator raising 
and lowering operators on general modified oscillator states. 

From the identities (|60H63() . we obtain the results: 

Claim 2: For A > 1 

rK x v {r) = VX + u-in^ir) + W+l 72£+i(r), (64) 
rft A (r) = V\~+^nl +1 {r)+^nlt\{r), (65) 

r f^o V M!T(A + zy) ^ 

OO 



^(D = g(-i)-'J ^;: ) 1) »rw. TO 
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(a) (b) 

FIG. 2: The action of the raising and lowering operators of radial wave functions compared to the actions of the standard 
raising and lowering operators of harmonic oscillator states. Multiplets of harmonic oscillator states of a given energy level are 
shown as horizontal lines in (b) labelled by radial v and SO(iV) angular momentum v. 



dr + v - 1 



-(^-^EHrtPS^M, (68) 



2 ^o V 
dr vA + v 



fl—V-\-l y \ f / 



(69) 



/'= 

The first of these equations follows from equations l|61|) and l|62[l and the observation that 

2r = A+(-A + |) +A(X- §) . (70) 

The second follows similarly from equations l|60(l and l|63|l with 

2r = A\\-\)+A{-\+\). (71) 

Equation l|66|) is obtained by expressing equation l|65|) as a recursion relation 



which is readily solved to give the desired result. Equation (I67[) is similarly obtained from equation H64|l . The last 
two equations of the claim are obtained from the identity 



A{X)-A^-X) = 2(± + X) 
\ar r J 



(73) 



Setting X = A — 3/2, and using equation i|66|l leads to equation (|68[1 . Setting X = —A + 1/2. and using equation 16 7|) 
leads to equation (|tj§l) . 



C. An 0(iV)-parity quantum number 

The above results show that radial matrix elements of type {Xfx\r 2 \Xv} , (Xn\l/r 2 \\v) , and (\fi\d 2 /dr 2 \\is) are 
obtained by use of the SU(1,1) Lie algebra and those of type (A± l,(j\r\\p), (A± 1, fi\l/r\\v), and (A± 1, ^i\d/ dr\Xv) 
are obtained by use of the factorization method. However, neither of the methods presented gives matrix elements of 
the type (\fi\r\\is), / r\\v) , and {\[i\d / dr\\v) . This is because odd powers of r and dj dr do not occur alone in the 
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space of polynomial functions of the basic {xi,pi} observables. For example, whereas r 2 is the quadratic r 2 = ^ x 2 , 
there is no polynomial expression for r = i/X^i x \- However, terms linear in r and d/dr do occur in combination with 
orbital functions. For example, when expressed in terms of M. N spherical polar coordinates, the Euclidean coordinates 
{xi} are of the form Xi = rQi(w), where Qi is proportional to a v = 1 SO(N) spherical harmonic. 

The implications of this observation are conveniently summarized in terms of an 0(iV)-parity quantum number 
7r = (— l) v associated with every SO(iV) irrcp of angular momentum v. It is seen that even and odd functions of the 
{xi} coordinates have even and odd 0(iV)-parity, respectively. It follows that the matrix elements (X' p,; i>'m|£i| Xv; vn) 
and (A'/j; v' m\pi\Xv; vn) vanish unless (— l) v = (— l) v+1 . This means that, in the evaluation of polynomial observables, 
the matrix elements of r and d/dr, for example, are only needed between states of opposite parity. Thus, the results 
of Sections IIV Al and IIV Bl lead to algebraic expressions for matrix elements of (positive and negative) integer powers 
of r and d/dr if one chooses basis states {\Xv; vm)} with A related to v by 

A, +1 = A„ ± 1. (74) 

The standard relationship A„ = v + N/2, given by equation (|17|) for the harmonic series of SU(1,1) irreps, obviously 
satisfies this condition. However, there are many other possibilities including, for example, having just two SU(1,1) 
irreps, one for even- and one for odd 0(7V)-parity states such that 

A+ = A_ ± 1 . (75) 

It should be ephasized that Hamiltonians that are not polynomials in {xi,pt} are also of interest. The Hamiltonian 
of the hydrogen atom, with a singular 1/r potential, is a standard example. However, it would appear that, for such 
Hamiltonians, the radial wave functions of the modified oscillator do not provide the most appropropriate basis wave 
functions. 



V. ORBITAL MATRIX ELEMENTS 

This section gives SO(iV)-reduced matrix elements of the basic v = 1, SO(iV) tensor, Q, defined by the expression 

x l = rQi (76) 

of the Cartesian coordinates in SO(iV) spherical polar coordinates. To make use of these reduced matrix elements one 
will, in general, need access to SO(iV) Clebsch-Gordan coefficients. The required CG ocefficients are known in SO (3) 
D SO(2) coupled bases for N < 6 (cf. Concluding remarks). 

A. SO(iV)-reduced matrix elements and a symmetry property 

SO(iV)-rcduced matrix elements are defined by the Wigner-Eckart theorem 

(v3™3|2m 2 K™-l) = (Vimi,lm 2 |u3"l 3 ) («3||Q||«l), (77) 

where {v\m\, 17712^37713) is an SO(iV) Clebsch-Gordan coefficient and \vm) is an orbital basis state whose wave function 
is an SO(iV) spherical harmonic y vm - It is then useful to define the SO(iV)-coupled action of a tensor operator, e.g. 
2, by 

[Q ® \vi)] V3ms = 2J \ v 3 m 3)( v 3 m 3\Qrn 2 \vimi)(vimi,lm 2 \v3m, 3 ) 

mi?ri2 

= \v 3 m 3 )(v 3 \\Q\\ Vl ). (78) 
An application of these definitions leads to the useful expression 

J [y„ 3 M ® C®y«,(w)]o du = k V3 (v 3 \\Q\\ Vl ), (79) 

and, hence, the identity 

k V3 (v3\\Q\\vi)=k Vl (v 1 \\Q\\v 3 ), (80) 
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where 

K = j [y v {uj)®y v {uj)] duj. (8i) 

The value of k v can be inferred, to within a phase factor, by making the expansion 

D^im = C nm y* m , (82) 

n 

which is always possible because the space £ 2 (Sn) has a real basis. Because spherical harmonics are defined to be an 
orthonormal basis of orbital wave functions, it follow that 

£|C„ m | 2 = l, (83) 

n 

and that 

k v = C nm ,{vm,vn\W) f y* m ,(u>)y vm (u) du = ^C nm (vm,vn\QQ). (84) 

m'ran mn 

Now, from the two identities 

Y \C nm \ 2 = d(v), vn\00) 2 = 1, (85) 

mn ran 

where d(v) is the dimension of the SO(N) irrep v, we conclude that 

Cnm = k v (vm, vn\ 00) and \k v \ 2 = d(v). (86) 

It follows that 



k v = J*Wy/d(y), (87) 

where cj>(v) is a phase angle, and 



(vs\\Q\\vx) = e i (*(»0-*(«.))^g) ( Vl \\Q\\v 3 ) . 

It is important to note, that reduced matrix elements are sometimes defined by expressing the Wigncr-Eckart 
theorem in the form 



/ in i \ t 1 I \ ^ z \\Q\\vx) 
{v 3 m 3 \Q m2 \vimi) = {virrii ,lm 2 \v 3 m 3 ) 



This adjusted definition of the reduced matrix elements results in a simplification of the symmetry relation i|88|) to 

(v 3 \\Q\\vi) = e *W*0-*(«.)) ( Vl \\Q\\v 3 ) . (89) 

B. Basic reduced matrix elements 

Reduced matrix elements of the basic v = 1, SO (AT) tensor, Q, are determined by the expansion 

x i = ^=(4 + c i ) (90) 

and the matrix elements of the harmonic oscillator raising and lowering operators. This expression shows that Q has 
non-zero reduced matrix elements (i>'||Q||i>) only for v' = v ± 1. 
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Let 

a t = _L (c t +ic t) (91) 

denote the highest weight component of the harmonic oscillator raising operators in which the operators {cj} of 
equation I|1U|) are regarded as Cartesian components of a v = 1, SO(N) tensor. Then, with A,„ = v + N/2, the 
harmonic oscillator ground state is the state |0) = |Ao^ = 0;v = m = 0) and a subset of excited states is given by 

\X v 0;vv) = -L^YlO), v = 0,l,2,--- (92) 
vvl 



The matrix element 

(A„+i0; v + 1, v + l[a. + | A„0; vv) = Vv+1 (93) 

then implies that 



(A, +1 0;« + l||rQ||A t ,0;w) = -j={\ v+1 0;v + l\\J\\\ v 0;v) = \J^~- (94) 

Thus, wc obtain the following result: 

Claim 3: SO(/V)-reduced matrix elements of the basic v = 1, SO(N) tensor, Q, defined by the expression Xi = rQi 
are given by 



WISH - ^ ^> +1 + ^<"-«-'» (95, 



The first term of this expression is obtained by factoring out the radial matrix clement (X v + 1, 0|r| A^O) = \f\, = 
y/v + N/2 in equation l|94|) . The second term is obtained from the first by use of the symmetry relationship (I88|) . 

For N = 3, for example, the customary expression of the Cartesian coordinates in terms of (r, 9, tp) spherical polar 
coordinates 

Xi = r sin# cos tp, 

X2 = r sin 9 sin tp, (96) 
£3 = r cos 9, 

defines the Cartesian components {Qi} of the TV = 3 tensor Q. The dimension of an SO(3) irrep of angular momentum 
I (a standard symbol for SO (3) angular momentum) is given by d(l) = 21 + 1 and, with the standard phase convention 
for SO(3) spherical harmonics, we obtain e 1 ^"' = (—!)'■ Thus, claim 3 gives 



o=a/^^ +1 +V^^- 1 - (97) 



For N = 5, the Wcyl dimension formula gives 



d(v) = ^{v + l){v + 2)(2v + 3) (98) 



and, with the phase convention e 1 ^^ = 1 (used in Ref. ^i|), claim 3 gives 



Wl+ vSi (99) 
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VI. COMBINED MATRIX ELEMENTS 



Reduced matrix elements of x = rQ are obtained by combining the expression for r, given in terms of radial raising 
and lowering operators by equations (|7L)|I and l|71|l . with those for the matrix elements of Q, given by claim 2. In 



particular, in the harmonic oscillator basis, for which A„ = v + N/2 (cf. fig.[21), one obtains: 

(\ v+1 v;v + l\\x\\\ v v;v) = ^(X v+1 v\A^ (X v - |)|A„^) (v + l\\Q\\v), (100) 

(\ v+1 ,v-l;v + l\\x\\\ v v;v) = ±(\ v+1 v\A(-\ v + |)|A„j/) (v + l\\Q\\v), (101) 

(A„_i, v + l;v- l\\x\\X v u; v) = ±(A„_i, u + l^-A. + |)|A,, v) (v - l\\Q\\v), (102) 

{\ v - lV -v- l\\x\\\ v u-v) = \{X v _ lV \A{X v - \)\K,v) («-l||e||w). (103) 

These identities have an immediate generalization. 



Claim 4: If a and (3 index any radial wave functions, the SO(iV) reduced matrix elements of the harmonic oscillator 
raising and lowering operators between states of SO (TV) angular momentum v and v' are given by 

(m/||c+||/3«) = ^ [{a\A*(v + |JV - \)\(3) S v ,, v+1 + {a\A*(-v - \N + |)|/3) 

x(v'\\Q\\v), (104) 
(av'\\c\\/3v) = ^ [(a\A(-v - |JV+ S v ', v+ i + (a\A(v + |JV - |)|/3) ^, v _ x ] 

x<t/||fi||«). (105) 



With the substitution Xi = -i=(c| + c^), these results follow immediately from equations H100HK)3|) when the radial 

wave functions are those of the harmonic oscillator basis. However, because the radial wave functions \V2^\v = 
0, 1,2,.. .} span the space of radial wave functions for any value of A, and because equations 1104|) and l|105(l with 
\/3v) = |A„i/; v) and \av') = \\ v iv'] v') hold for all values of v and v' they also hold for any radial wave functions. 
Because of the way the results of claim 4 are derived, they clearly incorporate the identity 



(av'\\x\\(3v) = (a\r\P) (v'\\Q\\v). 
What is more significant is that they also give the reduced matrix elements 



,v + l\\\p\\\Pv)=-ih(a\ 



(a, v — 1 



-dr 
d_ 
dr 



in- 

2 



r 

2 iv 



r 



(« + l||C||«> 
\y) («-l||ei|»> 



for the momentum operators 



p ^- lh dxrr-T2 {Cl 



4). 



(106) 

(107) 
(108) 

(109) 



VII. APPLICATION TO CENTRAL FORCE PROBLEMS 



The characteristic property of a central force problem is that its Hamiltonian is SO(./V)-invariant. For example, for 
a diatomic molecule the relevant Hilbert space is the space of £ 2 (R 3 ) wave functions in the relative coordinates of 
two atoms. For a free molecule in an isotropic space, the Hamiltonian is invariant under SO (3) rotations. Thus, its 
Schrodinger equation reduces to an equation in a single radial variable. Moreover, as already noted, the Hilbert spaces 
of radial wave functions arc independent of N. Thus, methods developed for the solution of central force problems in 
R 3 apply more generally. 

Consider, for example, the Hamiltonian of the quartic oscillator 



H= --V 2 + r 4 . 



(110) 
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Its spectrum was derived by Bell et al. 27] in two- and three-dimensional spaces by diagonalization in a basis of 
eigenstates of the harmonic oscillator Hamiltonian 

^HO = -^V 2 + iaV (111) 

with o=l. It has been considered more recently in five-dimensional space pfij . The addition of quartic terms, which 
lie in the SU(L1) enveloping algebra, to a nuclear collective model Hamiltonian has also been considered by several 
authors, e.g., [23. Some low-lying energy levels of the quartic oscillator in three dimensions are shown in Fig.|3 In 
repeating the calculations of ref. (231 for all L < 6 energy levels below 250 (in oscillator units), it was found that 21 s 
of computer time were needed for each L to compute these energy levels to an accuracy of 1 part in 10 12 when a = 1 
but only 2.6 s were needed with a = 1.6. 



20. 220,849,4 64,1 
L=0 



19. 915,973,0 18,2 

L=2 19. 217,580,4 82,7 
L=4 



18. 152,410, 722,0 
L=6 



10.099,944,419,6 

— L^j - 16. 046,193, 351,2 
L=3 



15. 081.646, 864,4 
L=5 



13. 379,336,5 52,6 
L=0 



13. 004,562,9 54,6 

L = 2 12.159,017,182,7 



L=4 



10. 099,944, 419,6 

Lli - 9. 401,160,1 55,8 

L=3 



7.335,729,995,2 

— rrrT- 6.830,307,934,3 

U L=2 



4. 478,039,2 19,52 
L=1 

2. 394,364,4 01 ,6 
L=0 



FIG. 3: Low energy levels of the quartic oscillator in three dimensions. 



A much larger gain is realized for the low-energy states of potentials with a minimum at a non-zero value of r as the 
following application to the study of a phase transition in the nuclear collective model shows. Such a calculation was 
performed originally j2|| in a five-dimensional harmonic oscillator basis. However, subsequent studies ^3] showed that 
it can be carried out much more efficiently using modified oscillator bases and the algebraic matrix elements associated 
with them. This gain in efficiency is particularly relevant, as discussed below, because it makes it practically possible 
to execute more sophisticated calculations in which the radial and orbital degrees of freedom are coupled. 

The Hilbert space for the nuclear collective model is £ 2 (M 5 ) and the Hamiltonian used in Ref. j2|| was 

H(a) = ~^V 2 + hi [(1 - 2a)r 2 + « r 4 ] , (112) 

where r is a radial coordinate for R 5 and M a mass parameter. This Hamiltonian is interesting because, as a passes 
through the critical value of 0.5, a phase transition occurs from a spherical vibrational phase, corresponding to a 
minimum value of V(r) at r — 0, to a rotational-vibrational phase for a > 0.5, corresponding to a minimum value of 
V(r) at yj (2a — l)/(2a) (in oscillator units proportional to \/A7). 



A. Variational calculations 



The benefits from the new develoments are maximized by selecting basis wave functions that are as close as possible 
to the eigenfunctions. For each value of a in the Hamiltonian, an optimal orthonormal basis of states \nvm) with 
wave functions of the form 



.(r.^V^M^P), (H3) 
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is obtained by varying the parameters a v and A„ for each value of the SO (TV) angular momentum v to minimize the 
energy expectation values of (0vm\H(a)\0vm) . The energy-levels obtained from the expectation values 

E nvm (a) w (nvm\H(a)\nvm), (114) 

with the u-dependent variationally-dctcrmincd values of a v and A„ , are shown as functions of a in comparison with 
accurately computed values in Fig. The close correspondence between the results and those computed numerically 
by diagonalization is remarkable. Even in the transition region, where restriction to single basis states is least 
successful, it does surprisingly well. 




(a) (b) 

FIG. 4: Comparison of the low-lying energy levels for the Hamiltonian H(a) with M = 100 as computed (a) by diagonalization 
and (b) by taking expectation values {nvm\H (a)\nvm) in single basis states with variationally chosen parameters as described 
in the text. The numbers shown on the right hand side of each figure are the values of the quantum numbers (vi, 1)2, • ■ • ; v) for 
each of the levels plotted. More detailed information is given for a — 1.5 in figureQ (Figure (a) was computed by P.S. Turner 

M-) 

Why the variational calculations with single basis wave functions are so good is illustrated in Fig. [S] which shows 
that a potential of the form 

V(r) = -br 2 +cr 4 , (115) 
which has a minimum at r = ro , is well fitted in the neighbourhood of its minimum by a potential 

W(r) =W + j^ + {ar) 2 , (116) 

with a scale parameter a chosen to give W(r) the curvature of the potential V(r) at ro. The radial eigenfunction for 
the ground state of a Hamiltonian with potential W(r) is given by T2.Q°(r) (in suitable units of r) and found to be an 
excellent approximation to that for the potential V(r). 

The values of a and A which minimize the variational ground-state energy are shown in Fig. |5| The figure shows 
that the wave function is essentially that of a spherical vibrator for a < 0.5 but that its width increases with a until 
the point at which the curvature of the potential V(r) = [(1 — 2a)r 2 + ar 4 ] at its r — minimum vanishes. 
With further increase in a, the wave function becomes that of a rotor-vibrator with equilibrium deformation given by 
r = [(A-l) 2 -9/4] 1 /4 /a 



B. Results from numerical diagonalization 

With the analytical expressions given in Sect. IIV Al variational calculations of the type shown above are easy to 
implement and give insightful information about the states of the system. Moreover, they are readily extended to 
obtain more precise results by diagonalization of the Hamiltonian in an optimal basis. 
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FIG. 5: Comparison of the potentials V(r) and W(r), as denned in the text, for fitted values of A and a. {TZq } is the radial 
wave function for the ground state of the Hamiltonian with potential W(r). 




iii i i 

0.5 1 1.5 2 



FIG. 6: Values of A and the inverse width parameter a/a(0), where a(0) = \/M, for the variational ground-state wave functions 
used in Fig. 0] (b) . 

If the objective is to obtain the energy eigenstates for the n lowest-energy states as accurately as possible when 
computed in a space spanned by the first N states of an ordered basis, then a simple prescription for selecting a near 
to optimal basis among a set of possible choices is to select the basis which minimizes the computed energy for the 
highest among the n lowest energy states. 

Other prescriptions can obviously be defined such as, for example, minimizing the sum of the energies of the lowest 
n states. However, it is not clear that much can be gained over the above simple prescription. 

For a > 1.5, for which the model nucleus described by the Hamiltonian H(a) of equation l|112|) is clearly non- 
spherical for M = 100, it was found that the energy of the lowest energy state of each SO(5) angular momentum 
v was obtained to within 1% accuracy with just one modified oscillator basis wave function. To achieve this level 
of accuracy with a conventional harmonic oscillator basis requires > 22 harmonic oscillator basis wave functions 
(depending on the value of a). Fig. [7| shows the excitation energies of the three lowest-energy states for each v < 6 
computed precisely and with just 5 basis states. 

The basis states were chosen for the v = states according to the above given prescription. This prescription could 
have been followed for each value of v. However, in practice, one will want to be able to compute electromagnetic 
transition matrix elements. Having fixed A for the v = states (which turned out to be 57) we then chose A = 57 for 
all even v and A = 58 for all odd v for the reasons discussed in Sect. IIV CI 

To obtain parallel results to the same level of accuracy in a harmonic oscillator basis requires ~ 25 basis states 
and a computation that takes ~ 10 times as long. For practical purposes, a reduction in the required number of 
radial basis states is particularly advantageous in treating systems in which there is a coupling between the orbital 
and radial degrees of freedom; this is because the product of the needed number of radial and orbital basis states can 
reach a large number. 
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(approx.) 



FIG. 7: Energy-level spectrum of the Hamiltonian H(a) for a = 1.5 and M = 100 for the three lowest-energy states of SO(5) 
angular momentum v < 6. The precise energies were computed by diagonalization with 100 basis states. The approximate 
values were computed with just 5 basis states chosen as described in the text. 



VIII. MORE GENERAL APPLICATIONS 

A. An N — 3 example 

In many problems of interest, the SO(N) rotational invariancc is broken. For example, in M 3 , it is of interest to 
study splittings of molecular energy levels by a crystal field. This might be done with a Hamiltonian of the type 

H=-^V 2 + V(r)+ X r 2 (3cos 2 9-l). (117) 

With the M = component of the L = 1 spherical tensor Q given, according to equation (|96|l . by Qio = cos 9, the 
function 3 cos 2 9 — 1 is proportional to the M = component of the L = 2 tensor (Q <g> Q) 2 , i.e., 

3 cos 2 (9- 1 = V6(Q® Q) 20 . (118) 
The reduced matrix elements of the tensor (Q (g> Q) 2 , obtained from those of Q by Racah recoupling, are given by 

(h\\(Q ® C)a||ii> = v/5(2Z + 1) W{h2, l 2 2; 1 2) (Z 2 || Q\\l) Q\\h). (119) 

It follows that 

(Z 2 m|(3cos 2 6>- l)\hm) =Y^m{2l + l){l 1 m 1 2Q\l 1 m)W{l 1 2M^l'2){h\\Q\\l){l\\Q\\h) 1 (120) 

where (hm, 20|Zim) is an SO(3) Clebsch-Gordan coefficient. These matrix elements are readily evaluated with the 
(l'\\Q\\l) matrix elements given by equation (|97|l . 

The effect on a spherical molecule of putting it into a quadrupole field is to deform it somewhat. However, if the 
molecule is already deformed the primary effect is to align the deformation of the molecule with the field. Thus, 
for a well-deformed diatomic molecule such as the HC1 molecule, an informative way to study the spectrum and 
cigenfunctions of the Hamiltonian i|117|) would be to start by solving for the Hamiltonian 

= -T^yV 2 + V(r) + X r 2 (3 cos 2 9-1). (121) 

This will give the splittings of energy levels due to the alignment effect. Solutions for the perturbed Hamiltonian 

H = H + X (r 2 - r 2 )(3 cos 2 0-1) (122) 

will then give the added perturbations of the spectra coming from the rotation-radial vibration coupling interactions 
and provide information about the rigidity of the molecule. 
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B. An N = 5 example 



For N > 3, there may be terms in the Hamiltonian which break the SO(N) invariance but retain rotational 
invariance with respect to a suitably defined SO(3) C SO(iV) subgroup. For example, a Hamiltonian of interest in 
the nuclear collective model is of the form 



1 



1 



H{a) = -—V 2 + -M[{l-2a) 



+ nr 3 {Q<g) Q® Q) Q , 



(123) 



where (Q® Q® Q)o is triple product of SO(5) D SO(3), v = 1,L = 2, Q tensors coupled to SO(3) angular momentum 
L = 0. This SO(3) coupled product is the L = component of a u = 3, SO(5) tensor. 
The matrix elements of (Q ® Q ® Q)o are given by 



E 



-1 



u i Q i L i 



(2Li + l)(2L 2 + l) 



W{L2L 2 L;L 1 2) 



2L+1 

x (u'a'i|| Q|jw2a2i2) (V2CX2L2W Q\\via\Li) {v\a\L\\\ Q\\vaL). 



(124) 



Note, however, that the reduced matrix elements appearing in this expression are SO(3)-reduced matrix elements; 
they are related to the SO(5)-reduced matrix elements of equation (j99(l by the expression 



{vxaxLxWQWvaL) = («aL,12||uiaiLi)(ui||Q||u), 



(125) 



where (vaL, 12|jwiaiLi) is an SO(5) Clebsch-Gordan coefficient in an SO(3) basis. An algorithm for computing such 
CG coefficients and tables of values has been given in Ref. [Tflj . 

The low-energy level spectrum of the Hamiltonian (|123J) . calculated with a large value of a jl^, is shown in Fig. 
The figure shows a ground-state rotational band and a sequence of excited (so-called gamma-) vibrational bands 
of the type given by the phenomenological Bohr-Mottelson nuclear collective model . 




FIG. 8: Low energy-level spectrum of the Hamiltonian of equation I)12HJ| with a large value of a and k = 50\/2/35. Each 
energy level is labelled by its SO(3) angular momentum; because of the interaction, the SO(5) angular momentum is no longer 
a good quantum number. Transition rates for electric quadrupole gamma-ray transitions indicated by arrows are shown in 
natural units beside the arrows; these transition rates are known in nuclear physics as reduced E2 transition rates. 
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IX. CONCLUDING REMARKS 



It has been shown that the matrix elements of polynomial Hamiltonians on a Euclidean space R can be determined 
algebraically to within SO (TV) Clcbsch-Gordan coefficients. 

The results obtained highlight the importance of algorithms for computing CG coefficients for various SO(iV) 
groups. Currently these coefficients are available for N < 6. Explicit expressions for SO(3) coefficients were given 
already in the 1931 edition of Wigner's book on Group Theory [3(| and can now be found in almost any book on 
angular momentum theory. The coefficients for SO (4), locally isomorphic to SO (3) x SO (3), are readily obtained 
from combinations of the SO (3) coefficients Q. An algorithm for calculating the needed coefficients for SO (5) has 
recently been given |19| in an SO(3) basis defined by regarding the fundamental 5-dimensional v = 1 irrep of SO(5) as 
carrying in irreducible L = 2 irrep of SO(3). This algorithm for SO(5) CG coefficients extends to SO(6) in an SO(5) 
D SO (3) basis by regarding the fundamental 6-dimensional irrep of SO (6) as spanning a sum of v = and v = 1 
SO (5) irreps and noting that the v = irrep of SO (5) is the trivial identity irrep. CG coefficients could undoubtedly 
be used with advantage for higher SO(iV) groups. For example, models for octupolc vibrations and rotations of nuclei 
or molecules with octupole deformations could be formulated on M 7 and would require SO(7) CG coefficients. 

It will interesting to see if the techniques introduced have a parallel extension to hydrogenic and modified hydrogenic 
systems in light of the recent developments of Fortunato and Vitturi [3lJ ■ 
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